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p Abstract 

A global next-to-leading order QCD analysis of unpolarized and polarized deep- 



inelastic scattering data is performed with parton distributions constructed in a sta- 
tistical physical picture of the nucleon. The chiral properties of QCD lead to strong 
relations between quarks and antiquarks distributions and the importance of the Pauli 
exclusion principle is also emphasized. We obtain a good description, in a broad range 
of x and Q 2 , of all measured structure functions in terms of very few free parameters. 
We stress the fact that at RHIC-BNL the ratio of the unpolarized cross sections for 
the production of W + and W~ in pp collisions, will directly probe the behavior of the 
d(x)/u(x) ratio for x > 0.2, a definite and important test for the statistical model. 
Finally, we give specific predictions for various helicity asymmetries for the W ± ,Z 
production in pp collisions at high energies, which will be measured with forthcoming 
experiments at RHIC-BNL and are sensitive tests of the statistical model for Au(x) 
and Ad(x). 



1 Introduction 



Deep-inelastic scattering (DIS) of leptons on hadrons has been extensively studied, 
over the last twenty years or so, both theoretically and experimentally. The principal 
goals of this physics program were, first to elucidate the internal proton structure, 
in terms of parton distributions, and more recently to test perturbative Quantum 
Chromodynamics (QCD), which generalizes the parton model. For the unpolarized 
structure functions, the advent of the HERA physics program gives us access to a 
broader kinematic range than fixed targets experiments, in x down to a few 10~ 5 
and in Q 2 up to several 10 4 GeV 2 , which allows testing perturbative QCD to next- 
to-leading order (NLO). As a result, the unpolarized light quarks (u,d) distributions 
are fairly well determined. Moreover, the data exhibit a clear evidence for a flavor- 
asymmetric light sea, i.e. d > u, which can be understood in terms of the Pauli 
exclusion principle, based on the fact that the proton contains two u quarks and 
only one d quark Q. Larger uncertainties still persist for the gluon (G) and the 
heavy quarks (s, c) distributions. From the more restricted amount of data on polar- 
ized structure functions, the corresponding polarized gluon and s quark distributions 
(AG, As) are badly constrained and we just begin to uncover a flavor asymmetry, 
for the corresponding polarized light sea, namely Au ^ Ad. Whereas the signs of 
the polarized light quarks distributions are essentially well established, Au > and 
Ad < 0, this is not the case for Au and Ad. The objective of this paper is to con- 
struct a complete set of polarized parton (all flavor quarks, antiquarks and gluon) 
distributions and, in particular, we will try to clarify this last point on the polarized 
light sea. 

The polarized parton distributions (PPD) of the nucleon have been extensively 
studied in the last few years [Q, |3| and in most models, the PPD are constructed from 
a set of unpolarized parton distributions, previously determined, from unpolarized 
DIS data. For example for each quark flavor qi(x), the corresponding Aq^x) is taken 
(at the input energy scale) such that 

Aqi(x) = di(x) ■ qi{x) , (1) 
2 



where cii(x) is a simple polynomial which has to be determined from the polarized 
DIS data. A similar procedure is used for antiquarks and gluons. As a result, the 
full determination of all unpolarized and polarized parton distributions involves a 
large number of free parameters, say around 20-25, which obviously shows a lack of 
simplicity. In addition, most of these models do not provide a flavor separation for 
the antiquarks q^x) and consequently for Aq^x). However, there are recent attempts 
to make this flavor separation, either using semi-inclusive polarized DIS data || or 
by means of a flavor-symmetry breaking ||. Our motivation for this work is to use 
the statistical approach to build up : qi, Aq iy q { , Aq { , G and AG, by means of a 
very small number of free parameters. A flavor separation for the unpolarized and 
polarized light sea is automatically achieved in a way dictated by our approach. 

The paper is organized as follows. In Section 2, we review the main points of our 
approach and we describe our method to determine the free parameters of the PPD 
with the set of experimental data we have used. In Section 3, we show the results 
obtained for the unpolarized DIS structure functions Ff ,d (x, Q 2 ) and xF^ N (x, Q 2 ) in 
a wide kinematic range, compared with the world data. We show the prediction of 
the ratio of unpolarized W + and W~ cross section at RHIC-BNL, which is sensitive 
to the d{x)/u(x) ratio, a challenging question for the statistical approach. Section 
4 is devoted to the polarized DIS structure functions gi ,d,n (x,Q 2 ). In Section 5, we 
give our predictions for single and double helicity asymmetries for the heavy gauge 
boson production (W , Z) in pp collisions at high energies, which are sensitive to Au 
and Ad and will be tested with forthcoming experiments at RHIC-BNL. We give our 
concluding remarks in Section 6. 

2 Basic procedure for the construction of the PPD 
in the statistical approach 

In the statistical approach the nucleon is viewed as a gas of massless partons (quarks, 
antiquarks, gluons) in equilibrium at a given temperature in a finite size volume. Like 
in our earlier works on the subject [|, [7|, ||], we propose to use a simple description of 
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the parton distributions p(x), at an input energy scale Q%, proportional to 



[exp[(x-X Qp )/x]±l]- 1 , (2) 

the plus sign for quarks and antiquarks, corresponds to a Fermi-Dirac distribution and 
the minus sign for gluons, corresponds to a Bose-Einstein distribution. Here Xq p is a 
constant which plays the role of the thermodynamical potential of the parton p and x 
is the universal temperature, which is the same for all partons. Since quarks carry a 
spin-1/2, it is natural to consider that the basic distributions are qf{x), corresponding 
to a quark of flavor i and helicity parallel or antiparallel to the nucleon helicity. This 
is the way we will proceed. Clearly one has = qf + q[ and Aqi = qf — q[ and 
similarly for antiquarks and gluons. 

We want to recall that the statistical model of the nucleon has been extensively 
studied in early and more recent papers in the literature || [Hj, but in these works at 



variance with our approach, the statistical picture is first considered in the nucleon 
rest frame, which is then boosted to the infinite-momentum frame. 

From the chiral structure of QCD, we have two important properties which allow 
to relate quark and antiquark distributions and to restrict the gluon distribution 
10: 



- The potential of a quark q\ of helicity h is opposite to the potential of the 
corresponding antiquark q~ h of helicity -h 

X 0q = ~ X 0q k ■ ( 3 ) 

- The potential of the gluon G is zero 

X 0G = . (4) 

From well established features of the u and d quark distributions extracted from 
DIS data, we anticipate some simple relations between the potentials: 

- u(x) dominates over d(x), therefore one can expect Xf u + Xq u > X$ d + Xq 8 

- Au(x) > 0, therefore Xq u > X$ u 

- Ad(x) < 0, therefore X^ d > X^ d . 



So we expect Xq u to be the largest thermodynamical potential and X^ d the small- 
est one. In fact, as we will see from the discussion below, we have the following 
ordering 

X 0u > X 0d ~ X 0u > x od ■ ( 5 ) 
Eq. (^|) is consistent with the previous determinations of the potentials [Q , including 
the one with dimensional values in the rest system [|TfJ. By using Eq. (3), this 
ordering leads immediately to some important consequences for antiquarks, namely 

i) d(x) > u(x), the flavor symmetry breaking which also follows from the Pauli 
exclusion principle, as recalled above. This was already confirmed by the violation of 
the Gottfried sum rule [[11], [L2|] . 

ii) Au(x) > and Ad(x) < 0, which remain to be checked and this will be done 
in hadronic collisions at RHIC-BNL (see Section 6). 

Note that since u + (x) ~ d + (x), we have 

Au(x) — Ad(x) ~ d(x) — u(x) , (6) 

so the flavor symmetry breaking is almost the same for unpolarized and polarized 
distributions. 

Let us now come back to the ordering in Eq. (5) to justify it. We consider the 
isovector contributions to the structure functions g\ and F2, which are the differences 
on proton and neutron targets. In the QCD parton model they read 

2xg { ^ n \x,Q 2 ) = \x[(Au + Au)(x,Q 2 ) - (Ad + Ad)(x,Q 2 )} ® AC N s(x,Q 2 ) , (7) 

and 

Ft n) (x,Q 2 ) = ± x [(u + u)(x,Q 2 )) - (d + d)(x,Q 2 ))]®C NS (x,Q 2 ) , (8) 

where ACns( x , Q 2 ) an d Cns{x, Q 2 ) denote the spin-dependent and spin- independent 
perturbative QCD coefficients 0. Since they differ only in a non- negligible way for 
very small 2, say x < 0.05, we see that 2xg ( f~ n ^ — i 7 2 P ~ n ' ) is only sensitive to the 
helicity minus components of the u and d quark distributions so we get 

[2xg { r n) -Ft n) ](x,Q 2 ) ~ -\[{u~ -d-)(x,Q 2 )) + (u- -T)(x,Q 2 ))]®C NS (x,Q 2 ) . 

(9) 
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Figure 1: The isovector structure functions 2xg[ p n \x) and n \x). Data is taken 
from Refs. H, U- 
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At this stage it is instructive to look at the available data shown in Fig. 1. We 
notice that these two functions have very similar shapes and their difference is small 
and mainly positive, except perhaps for large x. In order to try to identify the origin 
of this experimental fact, let us look at the integrals of these functions divided by x. 



The first one is twice the Bjorken sum rule ||13||, for which the best world estimate 



is Ibj = 0.176 ± 0.005 ± 0.007 and the second one is the Gottfried sum rule Jl I 



whose value is Iq = 0.235 ± 0.026 |]12| . As a result using Eq. (|9]) one obtains, say for 
Q2 = 5GeV 2 , 

f 1 dx[(d-(x) - u'ix)) + (T(x) - u~(x))} = 0.175 ± 0.06 . (10) 
J o 

Now from the NMC result on the Gottfried sum rule one has 

/ dx\d(x) - u(x)\ = 0.16 ± 0.03 . (11) 
Jo 

By comparing these two results we can assume, to a good approximation, the following 
relation for the helicity minus distributions 

d-(x) = u~{x) . (12) 

It follows from our procedure to construct antiquark distributions from quark distri- 
butions described above (see Eq. (3)), that we have automatically for the helicity 
plus antiquark distributions 

d + (x) = u + {x) , (13) 

which makes Eqs. (10) and (11) perfectly compatible. Indeed, as we will see below, 
Eq. (12) is rather well satisfied in the final determination of the distributions, after 
fitting the data. 

Let us now complete the description of our parametrization. As stated above, 
the essential ingredient for quarks and antiquarks is a Fermi-Dirac distribution, as 
shown in Eq. (2), but we expect this piece to die out in the small x region, so we 
have to multiply it by a factor AXQ q x b , where b > 0. In addition to A, a flavor and 
helicity independent normalization constant, we have introduced the factor Xq which 
is needed to get a good description of the data. It is not required by the simple Fermi- 
Dirac expression but, due to the ordering in Eq. (5), it will secure the correlation 
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between the shape of a given distribution and its first moment It is also in 

agreement with what has been found from data for the second and third moments of 
the valence partons fl6|] . The small x region is characterized by a rapid rise as x — > 
of the distribution, which should be dominated by a universal diffractive term, flavor 
and helicity independent, coming from the Pomeron universality. Therefore we must 
add a term of the form Ax b / [exp(x /x) + 1], where b < and A is a normalization 
constant. So for the light quarks q = u, dof helicity h — ±, at the input energy scale 
Ql = 4GeV 2 , we take 

AX h r b A^ h 

Qo) = F7 ^ / 1 + m , (14) 

y V '^ o; exp[(x -X^)/3c] + 1 exp(x/x) + l' V ; 



and similarly for the light antiquarks 

-hi _ ^(-^Og ) l%2b . (-, r\ 

Xq [X ' ^ ~ exp[(x + X - h )/x] + 1 + exp(x/x) + 1 ' 

Here we take 2b for the power of x and not b as for quarks, an assumption we will 
try to justify later. For the strange quarks and antiquarks, s and s, given our poor 
knowledge on both unpolarized and polarized distributions, we take the particular 
choice 

xs(x, Ql) = xs(x, Ql) = -[xu(x, Ql) + xd(x, Ql)] , (16) 

and 

xAs(x, Ql) = xAs(x, Ql) = \[xAd(x, Ql) - xAu(x, Ql)) . (17) 

This particular choice gives rise to a large negative As(x, Ql) and we will come back 
to it below, in the discussion of our results (see Section 4). The charm quarks c , 
both unpolarized and polarized, are set to zero at Ql = 4GeV 2 . Finally concerning 
the gluon distribution, as indicated above, we use a Bose-Einstein expression given 
by 

xG(x,Ql)= f G *\ , (18) 
exp(x/x) — 1 

with a vanishing potential and the same temperature x. This choice is consistent 
with the idea that hadrons, in the DIS regime, are black body cavities for the color 
fields. It is also reasonable to assume that for very small x, xG(x, Ql) has the same 



behavior as xg(x,Qo), so we wm take be = 1 + b. Since the normalization constant 
Aq is determined from the momentum sum rule, our gluon distribution has no free 
parameter. For the sake of completeness, we also need to specify the polarized gluon 
distribution and we take the particular choice 

xAG{x,Ql) = , (19) 

consistently with Eq. (4). As usual, the valence contributions are defined as q va i = 
q — q, so A and A are determined using the normalization of u va i(x) and d va i(x), whose 
first moments are respectively 2 and 1. 

To summarize our parametrization involves a total of eight free parameters 

x,X+ ,X " ,X M ,X+ ,6 ,6 and A . (20) 

In order to determine these parameters, we use a fitting procedure on a selection of 233 
data points at Q 2 values, as close as possible to our input energy scale Ql = 4GeV 2 
and the \ 2 value we obtain is 322. For unpolarized DIS, we have considered F% (x, Q 2 ) 
from NMC, BCDMS, E665 and ZEUS, F^(x,Q 2 ) from NMC, E665 and xF% N {x,Q 2 ) 



from CCFR [I2],fr7|-|p71. For polarized DIS we have considered gi ,d,n (x,Q 2 ) from 



SMC, E154 and E155 [14], ^9|. The five free parameters, one temperature and 
four potentials, which determine the Fermi-Dirac functions at the input energy scale 
Q 2 = 4GeV 2 are 

x = 0.09907, X+ = 0.46128, X ~ u = 0.29766, X ~ d = 0.30174 and X^ = 0.22775 . 

(21) 

From the above discussion which led us to Eqs. (12) and (13), we observe that the 
fit yields Xq u ~ X^ d , in agreement with our expectations. We show in Figs. 2, 3 the 
Fermi-Dirac functions Eq. (2) for light quarks and antiquarks, respectively. They 
exhibit a flat behavior in the small x region, at variance with the rising trend of the 
parton distributions, which ought to be described by the universal diffractive term 
(see Eqs. (14, 15)). It is also interesting to make one more observation from the above 
values of these potentials. It turns out that if we impose the following simple relations 
X^ d = Xq u , Xq u = 3/2Xq u and X^ d = 3/4X^5 we can get an equally good fit of the 
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data with Xq u = 0.30549. This choice reduces the number of free parameters for the 
potentials from four to one, but we cannot justify it and it might be fortuitous. 
For the remaining three free parameters b , b and A, the fit gives the values 

b = 0.40962, b = -0.25347 and A = 0.08318 . (22) 

Finally the parameters A, A and A G determined by normalization conditions and 
momentum sum rule, have the following values 

A = 1.74938, A = 1.90801 and A G = 14.27535 . (23) 

We show in Figs. 4, 5, the different helicity components of the light quarks and 
antiquarks, respectively, at Q 2 = 20GeV 2 , after a NLO evolution. They all have the 
same rising behavior in the small x-region, which is driven by the universal diffractive 
term. We note that after NLO evolution, we still have u~ ~ d~ and u + ~ d + . In Fig. 
6, we display the x-shapes of the full set of unpolarized parton distributions, where 
one sees a non zero c quark distribution generated by the Q 2 evolution. Note that 
u~ ~ d~ implies 

u(x) - d(x) ~ Au(x) - Ad(x) (24) 

at least up to Q 2 ~ 20GeV 2 . This is a specific feature of our quark distributions which 
is not fulfilled in most parametrizations @. In order to compute the evolved distri- 
butions, we have used the following method: all parton distributions and splitting 
functions are decomposed in terms of Chebyshev polynomials, then the DGLAP equa- 
tions at NLO are solved by a Runge-Kutta method |3(J or by a semianalytic method 



il| (the renormalization scheme adopted is MS with A.j^g[nf = 3] = 300MeV). We 
have checked that both methods are consistent, within the numerical accuracy. All 
the results shown in the figures are calculated at NLO. 

Let us now comment on the values obtained for some of these parameters. First, 
it is interesting to note that A and A, the normalizations of quarks and antiquarks, 
come close to each other. Next, concerning the power of x for which we took, b for 
quarks and 2b for antiquarks, let us try to understand this fact and the value of b. 
The universal diffractive term, involving b and A, is absent in the physical quantity 
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xd(x) — xu(x), which was extracted recently in Ref. [0 (see Section 3 for further 
comments on these data). This quantity rises in the region 0.01 < x < 0.2, whereas 
the antiquark Fermi-Dirac functions decrease (see Fig. 3). This rising behavior is 
simply obtained if one multiplies these functions by x ' 8 . Another physical quantity 
which does not contain the universal diffractive term is xF% (x), because it involves 
the difference qi(x) — <?j(x). In the region 0.01 < x < 0.2, quarks dominate over 
antiquarks (see Fig. 6) and we know that xF^ N (x) rises. This behavior is obtained 
if one multiplies the quark Fermi-Dirac function by x ' 4 , in agreement with Eq. (22). 

In the next section we will compare our calculations obtained using these parton 
distributions, with the existing experimental world data, for the unpolarized DIS 
structure functions F^' (x, Q 2 ) and xF£ N (x, Q 2 ), in a wide kinematic range. We will 
also show the prediction of the ratio of unpolarized W + and W~ cross sections at 
RHIC-BNL, which is sensitive to the d(x) /u(x) ratio, a challenging question for the 
statistical approach. 
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Figure 2: The Fermi-Dirac functions for quarks = XQ q /(exp[(x — Xq^/x] + 1) at 
the input energy scale Ql = 4GeV 2 , as a function of x. 
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Figure 3: The Fermi-Dirac functions for antiquarks = 1 / X^(exp[(x + X^) / x] + 1) 
at the input energy scale Ql = 4GeV 2 , as a function of x. 
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Figure 5: The different helicity components of the light antiquark distributions after 
NLO evolution, at Q 2 = 20GeV 2 , as a function of x. 
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3 Experimental tests for unpolarized DIS 



We first consider fj,p and ep DIS for which several experiments have yielded a large 
number of data points on the structure function F$(x, Q 2 ). We have compared our 
calculations with fixed target measurements NMC, BCDMS and E665, which cover 
a rather limited kinematic region in Q 2 and also with the data at HERA from the 
HI and ZEUS Collaborations. These last data cover a very large Q 2 range, up to 
Q 2 = 10 4 GeV 2 or so and probe the very low x region which is dominated by the 
rising behavior of the universal diffractive term. 

The comparison of our results with the data is shown in Figs. 7, 8, 9 and 10. We 
notice in Fig. 7 that the very low Q 2 range accessible by E665, requires a downwards 
Q 2 evolution which is achieved successfully down to Q 2 = lGeV 2 or so, but we fail 
to reproduce the data much below than that. The HI and ZEUS data are fairly well 
described as shown in Figs. 8 and 9, respectively. Finally, we present a compilation 
of the data, including NMC and BCDMS, in Fig. 10 which is in very good agreement 
with our theoretical curves. From measurements over a large Q 2 range, it is possible 
to improve the determination of the gluon density by analyzing the scaling violations. 
On Fig. 11 we see that xG(x, Q 2 ) exhibits a fast rising behavior in the low x region, 
which is fairly consistent with our simple parametrization (see Eq. (18)). Note that, 
xG(x, Q 2 ) is predicted to increase with Q 2 as shown in the figure. Next we consider 
F?f(x, Q 2 ) obtained on a deuterium fixed target from NMC, BCDMS and E665 data. 
The comparison of these very accurate data with our results is shown on Figs. 12, 13 
and 14, respectively. The agreement is also excellent, except for the very low Q 2 
region of E665, as for the proton case. Finally the high statistics vN DIS data from 
CCFR allows to extract the xFg N (x, Q 2 ) structure function, which is successfully 
compared to our results on Fig. 15. 

To complete our tests of the unpolarized parton distributions, we must come 
back to the important question of the flavor asymmetry of the light antiquarks. Our 
determination of u(x, Q 2 ) and d(x, Q 2 ) is perfectly consistent with the violation of the 
Gottfried sum rule, for which we found Iq = 0.2493 for Q 2 = 4GeV 2 . Nevertheless 
there remains an open problem with the x distribution of the ratio d/u for x > 0.2. 
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According to the Pauli principle this ratio should be above 1 for any value of x. 



However, recently the E866/NuSea Collaboration [32] has released the final results 
corresponding to the analysis of their full data set of Drell-Yan yields from an 800 
GeV/c proton beam on hydrogen and deuterium targets and they obtain the ratio, 
for Q 2 = 54GeV 2 , d/u shown in Fig. 16. Although the errors are rather large in the 
high x region, the statistical approach disagrees with the trend of the data. Clearly 
by increasing the number of free parameters, it is possible to build up a scenario 
which leads to the drop off of this ratio for x > 0.2. For example this was achieved 
in Ref. [3B|, as shown by the dashed curve in Fig. 16. There is no such freedom in 
the statistical approach, since quark and antiquark distributions are strongly related. 
One way to clarify the situation is, either to improve the statistical accuracy on 
the Drell-Yan yields which seems rather unlikely, or to call for the measurement of 
another observable sensitive to u(x) and d(x). One possibility is the ratio of the 
unpolarized cross sections for the production of W + and W~ in pp collisions, which 
will directly probe the behavior of the d(x) /u(x) ratio. Let us recall that if we denote 
Rw{y) — {da w+ j dy) / (da w /dy), where y is the W rapidity, we have |37|] at the 
lowest order 

7~> / M 2 1 = u(x a , M^)d(x b , Agj + djxg, M^)u(x b , Mgj 

nw[y ' w) d(x a ,]^)v(x b ,]^)+V(x a ,M&)d(x b ,M&)' { } 

where x a = \fre y ', x b = \fre~ y and r = M^/s. This ratio Rw, such that Rwiy) = 

Rw(—y), is accessible with a good precision at RHIC-BNL []38| and at ^fs = 500GeV 

for y — 0, we have x a = x b = 0.16. So R w (0, M^) probes the d(x)/u(x) ratio at 

x = 0.16. Much above this x value, the accuracy of Ref. |32|] becomes poor. In Fig. 

17 we compare the results for Rw using two different calculations. In both cases we 

take the u and d quark distributions obtained from the present analysis, but first 

we use the u and d distributions of the statistical approach (solid curve in Fig. 16) 

and second the u and d from Ref. [[3(J (dashed curve in Fig. 16). Notice that the 

energy scale is much higher than in the E866/NuSea data, so one has to take into 

account the Q 2 evolution. At a/s = 200GeV for y — 0, we have x a = x b = 0.40 and, 

although the W ± yield is smaller at this energy, the effect on Rw(®, ^w) is strongly 

enhanced, as seen in Fig. 17. This excellent test, must be done in the near future. 
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Figure 7: F%(x, Q 2 ) as function of Q 2 for fixed x, E665 data |^TJ. The function 
c(xi) = 0.6(19 — i), i = 1 corresponds to x — 8.9 10~ 4 . 
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Figure 8: F$(x, Q 2 ) as function of Q 2 for fixed x, HI data p3| , |33|. The function 
c(xi) = 0.6(19 — i), i = 1 corresponds to x = 1.78 10~ 4 . 
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Figure 9: F^x^Q 2 ) as function of Q 2 for fixed x, ZEUS data |24| , |25|j . The function 
c(xi) = 0.6(19 — z), i = 1 corresponds to x = 6.3 10~ 5 . 
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Figure 10: F%{x, Q 2 ) as function of Q 2 for fixed x, c(x) = O.Q(i x — 0.4), i x — 1 - 
0.32, rebinned data HI, ZEUS, E665, NMC, BCDMS. (Presentation of data, co 
of R. Voss). 
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Figure 11: Comparison of xG(x,Q 2 ) at Q 2 = 20 — 30GeV 2 (dashed-solid) with ex- 
perimental determination from NMC [lTj, HI |35j and ZEUS 0j experiments. 
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Figure 13: F^(x,Q 2 ) as function of Q 2 for fixed x, BCDMS data 0. The function 
c(xi) = 0.6(19 — i), i — 1 corresponds to x = 7 10~ 2 . 
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Figure 14: F$(x,Q 2 ) as function of Q 2 for fixed x, E665 data [21]. The function 
c(xi) = 0.6(19 — i), % — 1 corresponds to x = 8.9 10~ 4 . 
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Figure 15: xFg N (x, Q 2 ) as function of Q 2 for fixed x, CCFR data [27]. The function 
c(xi) = 0.6(19 — i), i — 1 corresponds to x — 7.5 1CT 3 . 
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Figure 16: Comparison of the data on d/u(x, Q 2 ) from E866/NuSea at Q 2 = 54GeV 2 
||32|| , with the prediction of the statistical model (solid curve) and the set 1 of the 
parametrization proposed in Ref. |3(| (dashed curve). 
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Figure 17: Theoretical calculations for the ratio Rw(y, ^w) versus the W rapidity, 
at two RHIC-BNL energies. Solid curve (y/s = 500GeV) and dashed curve (y/s = 
200GeV) are the statistical model predictions. Dotted curve (y/s = 500GeV) and 
dashed-dotted curve (y/s = 200GeV) are the predictions obtained using the d(x)/u(x) 
ratio from Ref. 
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4 Experimental tests for polarized DIS 



Since our approach is based on the direct construction of the quark and antiquark 
distributions of a given helicity qf and qf , from the previous results we immediately 
obtained Aqi and Aq { for each flavor. We display in Fig. 18 these distributions 
xAf(x,Q 2 ) versus x, at Q 2 = 20GeV 2 , after a NLO evolution. As we mentioned 
earlier, we took AG = As = As = at the input energy scale, but as shown in the 
Figure, this is no longer true after the Q 2 evolution. We notice that the distributions 
are positive for u, u and G and negative for d, d, s(s) and c, which remain extremely 
small. We have also checked that our Aq,i(x) satisfy the positivity conditions at the 
leading twist level obtained in Ref. [pSR . Recently, the HERMES Collaboration has 
presented new semi-inclusive data |^0] with greater precision, which allows a good 



flavor separation of the light quarks and these data are shown in Fig. 19. However 
from this data analysis one cannot achieve a flavor separation for the antiquarks, 
so Aq s /q s represents the polarization of sea quarks assuming flavor symmetry, i.e. 
Aq s /q s = Au/u = Ad/d = As/s. Our model calculations give a very good description 
of the u and d quark polarizations as shown in the Figure. The comparison between 
our three curves Au/u, Ad/d and As/s and the poor accuracy sea quark data does 
not allow one to draw any conclusion. 

In Fig. 20 we show a data compilation of the polarized structure functions 
gi' d,n (x, Q 2 ) from different current experiments on proton, deuterium and helium 
targets, evolved at a fixed value Q 2 = 5GeV 2 . The x dependence is in fair agreement 
with our results and we predict, in the small x region, a fast rising behavior for g\ and 
a fast decreasing behavior for g™, due to the antiquark contributions. This cannot be 
tested so far, due to the lack of precise data. The Q 2 dependence for fixed x values 
is displayed in Figs. 21 and 22. 

Here we would like to comment on the choice we made for As(x, Q 2 ) (As) in 
Eq. (17). Clearly this polarized quark distribution is very badly known and we have 
constrained its first moment by assuming the validity of the second Bjorken sum rule, 
namely 

Aq 8 = Au + An + Ad + Ad - 2(As + As) = 3F - D , (26) 
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where F and D are the hyperon beta decay constants, so that 3F — D = 0.579. At 
Q 2 = 4GeV 2 , we have found for the following first moments, Au = 0.6678, Au = 
0.0464, Ad = -0.2576, Ad = -0.0865 and As = As = -0.0443, so Ag 8 = 0.547. 
Notice that the first Bjorken sum rule is also very well satisfied, since we get I^j = 
0.1764. It is interesting to realize that the contribution of the antiquarks to Ibj, 
whose value in our case is 0.0221, is positive and relatively large because Am and Ad 
have the same signs as Am and Ad, respectively. 

Finally we turn to the important issue of the large x behavior of the polarized 
quark distributions. This kinematic region has been poorly explored experimentally 
so far, but there are different theoretical scenarios [f|lj], when x is near 1, for the 
asymmetries A v { ' n (x, Q 2 ), measured in polarized DIS. We recall the definition of the 
asymmetry A\{x, Q 2 ), namely 

A( n2 ,_ 9i(x,Q 2 ) 2x[l + R(x,Q 2 )} 

Al{X ' Q } ~ F 2 (x,Q 2 ) [l+7 2 (x,Q 2 )] ' { ° 

where 7 2 (x, Q 2 ) = 4M 2 x 2 /Q 2 and R(x, Q 2 ) is the ratio between the longitudinal and 
transverse photoabsorption cross sections. In the case where the u quark dominates, 
we get 

Au(x,Q 2 ) [l + R(x,Q 2 )] 
X ~ u(x,Q 2 ) [l + 7 2 (x,Q 2 )] ' 1 J 

When x — > 1 for Q 2 = 4GeV 2 , R is of the order of 0.30 or less and 7 2 (x, Q 2 ) is close 

to 1, so A\ ~ 0.6Au(x)/u(x). It seems unlikely to find A\ —>■ 1, unless one lets Q 2 

go to infinity or one violates positivity. We show in Fig. 23 a compilation of the 

world data for Ai' n (x, Q 2 ) at Q 2 = 4GeV 2 , with the results of our calculations up to 

x = 1, where we have Au/u = 0.77 , Ad/d = —0.46 and indeed we find A v { n < 1. 

This specific prediction should be confronted with the very accurate data on A\ in 



the large x region which is expected soon from Jefferson Lab |42 
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Figure 18: The different polarized parton distributions (/ —u, d, u, d, s(s), c and G) 
after NLO evolution, at Q 2 = 20GeV 2 , as a function of x. 
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Figure 19: HERMES data on (Am + Au)/(u + u), (Ad + Ad)/(d + d), Aq s /q s as 
function of x at fixed Q 2 = 2.5GeV 2 gOj. The curves are our model calculations. For 
the sea quarks Au/u (solid curve) , Ad/d (dashed curve) and As/s (dotted curve). 
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Figure 20: $ An (x, Q 2 ) as function of x for different Q 2 values, from E155, E154, 
E143, SMC, HERMES experiments. The curves correspond to our model predictions 
at Q 2 = 5GeV 2 . 
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Figure 21: 2 n g\(x ) Q 2 ) as function of Q 2 for different x values, n = corresponds to 
x = 0.75 and n = 16 to x = 7.5 1CT 3 . Experimental data are rebined to the nearest 
x values. 
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Figure 22: g™(x,Q 2 ) as function of Q 2 for different x values. The function c(xj) = 
19 — i, % — corresponds to x — 7.5 1CT 3 . Experimental data are rebined to the 
nearest x values. 
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Figure 23: Compilation of the asymmetries A\ and A" from E155, E154, E142, E143, 
EMC, SMC and HERMES experiments [43]- [45]. The curves correspond to our model 
predictions at Q 2 = 4GeV 2 . 
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5 Helicity asymmetries for weak boson production 
at RHIC-BNL 



Next we propose some tests of our PPD in hadronic collisions in the framework of 



the spin program at RHIC-BNL f38fl . As we have seen in Section 3, the production 
of in pp collisions is very relevant to probe the behavior of the d(x) /u(x) ratio 
and since W bosons are produced through a pure V — A interaction, they are also an 
ideal tool to study the spin structure of the nucleon. 

Let us consider the parity- violating helicity asymmetry A PV (W) defined as 

< V (W) = *?l*-*fl!% • (»> 



where ± stands for the helicity of one polarized proton beam. For W + , at the lowest 
order of the Drell-Yan production mechanism, it reads 



a pv (w+) = A^q, M^)d(x b , Mgj - Ad(x a , M^)u(x b , M|) 
L 1 '~ u(x a ,M^)d(x b ,M^) + d(x a ,M^)u(x b ,M^) ' 1 } 

assuming the proton a is polarized. Here defined as in Eq. ( p5[) and for W 

production one interchanges the quark flavor u and d. The calculation of these asym- 
metries is therefore very simple and the results using our PPD are presented in Fig. 
24 at yfs = 350GeV and 500GeV. The asymmetries decrease for increasing energy 
and we recall that higher-order corrections have very small effects on these predic- 
tions [47|. The general trend of A^ V (W) can be easily understood and, for example at 
^ = 500GeV near y = +1, A£ V (W + ) ~ Au/u and A P L V (W~) ~ Ad/d, evaluated at 

x = 0.435. Similarly for near y = -1, A PV (W + ) Ad/d and A P L V {W-) Au/u, 

evaluated at x — 0.059. Given the expected rates for W ± production at RHIC-BNL 
and the high degree of the proton beam polarization []38[| , it will be possible to check 
these predictions of the statistical approach to a high accuracy, in particular for the 
flavor separation of the antiquarks polarized distributions. It will provide the first 
reliable determination of Am and Ad. One can also consider the asymmetry A PV 
for the Z° production whose expression is given in Ref. [|B|. Our prediction for this 
asymmetry is displayed in Fig. 25, but since the Z° is not a pure left-handed object 
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the interpretation of the result in terms of the PPD is less obvious than in the 
case. Moreover the Z° production rate at RHIC-BNL will be less copious, so the 
expected precision will be reduced. 

In pp collisions for production where both protons beams are polarized, there 
is another observable which is sensitive to the antiquark polarized distributions, that 
is the parity-conserving double helicity asymmetry A^(W) defined as 



apc (w \ _ da ++/ d V + da^_jdy - dgW_/dy - da^ + /dy , o , 

LL[ } daf + /dy + daW_/dy + day_/dy + daW + /dy ' ' ' 



For W + production it reads, at the lowest order, 



{ pc (w +) = Au ( x *> M^)Ad(x b , M|j + Ad(x a , M^)Au(x b} 
LL[ ' ~ ~ u{x a , M^)d{x b , M&) + d(x a , M^)u(x b , Mfa) 



We simply notice that if the antiquarks are unpolarized, i.e. Aq^x) = 0, it leads 
immediately to A^(W + ) = 0. Clearly A^(W + ) is symmetric around y = and 
one obtains All(W~) by interchanging the quark flavor u and d. Given the signs of 
the PPD we have obtained in the statistical approach, it is obvious that A^iW^) 
are positive in the y range we are considering, as shown in Fig. 26. The asymmetries 
decrease for increasing energy but they remain sizeable even at the highest RHIC-BNL 
energy. 

We have also calculated this asymmetry for Z° production and the results are 
presented in Fig. 27. The sign change with respect to Af£'(W /± ) is due to the fact 
that the expression of |37| is driven by terms of the type — Au(x a ) Au(x b ) 



or — Ad(x a ) Ad(x b ), which are both negative in the statistical approach. 

Finally let us mention briefly a realistic possibility of having polarized 3 He at 
RHIC-BNL, which allows to consider pn collisions with polarized neutrons. The same 
asymmetries can be calculated for this case and, in particular, we have checked that 
A^ v (W^) for pit collisions is directly obtained from A^fW 111 ) for p~p* collisions, by 
exchanging W + and W~, as a consequence of isospin symmetry. For ~p*n collisions 
with only the proton beam polarized, the results are very close to those obtained for 
~]?p collisions. 
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Figure 25: The parity violating asymmetry A^ v for pp — > Z° production versus the 
Z° rapidity at y/s = 350GeV (dashed curve) and y/s = 500GeV (solid curve). 
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Figure 26: Parity conserving double helicity asymmetry for pp — > W ± production 
versus the W rapidity at = 350GeV (dashed curve) and y/s = 500GeV (solid 
curve). 
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Figure 27: Parity conserving double helicity asymmetry A^l for pp — > Z° production 
versus the Z° rapidity at y/s = 350GeV (dashed curve) and y/s = 500GeV (solid 
curve). 



43 



6 Concluding remarks 



We have constructed a new set of polarized parton distributions in the framework 
of a statistical physical picture of the nucleon. The chiral properties of QCD lead 
to simple relations between quark and antiquark distributions. We have obtained 
a good description of all unpolarized and polarized structure functions measured in 
DIS, Fj ' (x, Q 2 ), xFg N (x, Q 2 ) and gl ,d,n (x, Q 2 ), in a wide kinematic range, in terms of 
eight free parameters, a small number which could be reduced even further. We have 
also proposed a simple expression for the gluon distribution G(x, Ql), at the input 
scale Qq, with no additional free parameter. In view of the poor present knowledge 
of the polarized gluon distribution we took AG(x, Ql) = 0. Our analysis of the data 
was done by using a DGLAP Q 2 evolution at NLO, to get access to a broad kinematic 
range of x and Q 2 . This approach predicts specific flavor symmetry breaking of the 
light antiquarks. First the inequality d(x,Q 2 ) > u(x,Q 2 ), a natural consequence of 
the Pauli principle, which is questioned for x > 0.2 by the E866 experiment |]32|j . 
To answer this challenging issue, we propose to measure the ratio of the unpolarized 
cross sections for the production of W + and W~ in pp collisions at RHIC-BNL, for 
which we have a definite prediction. Second, this approach leads to Au > and 
Ad < 0, so these first moments of the light antiquarks distributions give a positive 
contribution to the Bjorken sum rule. The precise tests for the x and Q 2 dependence 
of these distributions will also be done at RHIC-BNL, by the measurements of the 
helicity asymmetries in the W ± and Z° production, which we have calculated. The 
strange quarks unpolarized and polarized distributions, are constructed by means of 
empirical expressions in terms of the light antiquark distributions, but this should 
be improved in the future with some more fundamental arguments. We have also 
stressed the relevance of the large x region for the light quark distributions and the 
behavior of our distributions can be tested in current experiments at Jefferson Lab. 
Although we still miss a deep understanding of all the features of this approach, as 
well as a clear interpretation of some of the outgoing parameters, we think it is very 
promising and challenging in view of the future tests we have identified. 
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